Anisotropic in-plane optical conductivity in detwinned Ba(Fei_2;Co2;)2As2 
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We study the anisotropic in-plane optical conductivity of detwinned Ba(Fei-xCox)2As2 single 
crystals for x=0, 2.5% and 4.5% in a broad energy range (3 meV-5 eV) across their structural and 
magnetic transitions. For temperatures below the Neel transition, the topology of the reconstructed 
Fermi surface, combined with the distinct behavior of the scattering rates, determines the anisotropy 
of the low frequency optical response. For the itinerant charge carriers, we are able to disentangle 
the evolution of the Drude weights and scattering rates and to observe their enhancement along the 
orthorhombic antiferromagnetic a-axis with respect to the ferromagnetic 6-axis. For temperatures 
above Ts , uniaxial stress leads to a finite in-plane anisotropy. The anisotropy of the optical conduc- 
tivity, leading to a significant dichroism, extends to high frequencies in the mid- and near-infrared 
regions. The temperature dependence of the dichroism at all dopings scales with the anisotropy 
ratio of the dc conductivity, suggesting the electronic nature of the structural transition. Our find- 
ings bear testimony to a large nematic susceptibility that couples very effectively to the uniaxial 
lattice strain. In order to clarify the subtle interplay of magnetism and Fermi surface topology we 
compare our results with theoretical calculations obtained from density functional theory within the 
full-potential linear augmented plane- wave method. 

PACS numbers: 74.70.Xa,78.20.-e 



I. INTRODUCTION 

In many unconventional superconductors including 
cuprates and Fe-based pnictides, superconductivity 
emerges from a complicated soup of competing phases in 
the normal state when magnetism is suppressed by dop- 
ing, pressure or other external parameters [U |2] . This 
multiplicity of phases includes nematicity, defined as the 
spontaneously broken C4 rotational symmetry of the 
square lattice, and a novel form of magnetism, arising 
from either orbital currents or antiferromagnetic fluctu- 
ations. In the Fe-based pnictide superconductors (for a 



review with a comprehensive references list see Ref. [3]), 
nematic correlations and antiferromagnetic fluctuations 
have also been recently connected with symmetry break- 
ing competing phases. Inelastic neutron scattering ex- 
periments revealed anisotropic magnetic excitations [4], 
coupled to the structural tetragonal-orthorhombic transi- 
tion at Ts . This structural transition elongates the Fe-Fe 
distance in the ah plane along the a-axis direction and 
contracts it along the perpendicular 6-axis. It turns out 
that below the Neel temperature T/v (< Ts) the spins 
present ferromagnetic correlations along the shorter or- 



thorhombic 6-axis and antiferromagnetic stripes along 
the longer a-axis [4 . The two- fold anisotropy is also 
evident in quasiparticle interference observed via scan- 
ning tunnelling microscopy measurements [5 and further 
confirmed by angle-resolved-photoemission-spectroscopy 
(ARPES) data, collected on crystals for which the in- 
cident beam size was comparable to the size of a single 
structural domain [6 . Furthermore, quantum oscillations 
in the parent compound reveal that the reconstructed 
Fermi-surface (FS) comprises several small pockets [7]. 
The smallest of these pockets is essentially isotropic in 
the a6-plane, but the other, larger pockets are much more 
anisotropic. 

The first ARPES observations of an anisotropic elec- 
tronic dispersion [6 motivated an intensive research ac- 
tivity also with probes for which any impact of the elec- 
tronic anisotropy would be obscured by the formation of 
dense twin domains. These correspond to adjacent mi- 
croscopic domains as small as a few microns with alter- 
nating orthorhombic a and b axes [8 . Two distinct meth- 
ods have been employed so far to detwin the specimen: 
application of uniaxial stress |8J |9J and of an in-plane 
magnetic field [10]. The former method, employed here, 
is superior in order to achieve an almost complete detwin- 
ning. Recent ARPES measurements [1JJ112J of detwinned 
single crystals of Ba(Fei_2;Coa;)2As2 reveal an increase 
(decrease) in the binding energy of bands with dominant 
dyz {dxz) character on cooling through Tg [11 , leading to 
a difference in orbital occupancy. The splitting of the dxz 
and dyz bands is progressively diminished with Co sub- 
stitution in Ba(Fei_a;Coa;)2As2, reflecting the monotonic 
decrease in the lattice orthorhombicity [2{a — b)/{a-\-b)]. 
For temperatures above Tg , the band-splitting can be in- 
duced up to rather high temperatures by uniaxial stress. 

Mechanically detwinned crystals also provide a suit- 
able playground in order to explore the intrinsic in-plane 



anisotropy of the transport properties. Measurements of 
the dc resistivity as a function of temperature of the sin- 
gle domain parent compounds BaFe2As2, SrFe2As2 and 
CaFe2As2 (i.e., so called 122 iron pnictides) reveal a mod- 
est in-plane dc anisotropy for temperatures below Tg, 
with the resistivity in the ferromagnetic direction larger 
than along the antiferromagnetic direction [9l |T3l IT4] . 
Substitution of Co, Ni or Cu suppresses the lattice or- 
thorhombicity [15 , but in contrast the in-plane resis- 
tivity anisotropy is found to initially increase with the 
concentration of the substituent, before reverting to an 
isotropic in-plane conductivity once the structural transi- 
tion is completely suppressed [H [16] . Perhaps coinciden- 
tally, the onset of the large in-plane anisotropy for the 
cases of Co and Ni substitution occurs rather abruptly 
at a composition close to the start of the superconduct- 
ing dome. 

For temperatures above T^, there is a remarkably large 
sensitivity to uniaxial pressure, leading to a large in- 
duced in-plane resistivity anisotropy that is not observed 
for overdoped compositions [9 . There is no evidence in 
thermodynamic or transport measurements for an ad- 
ditional phase transition above Tg for unstressed crys- 
tals, implying that the induced anisotropy is the result 
of a large nematic susceptibility, rather than the pres- 
ence of static nematic order. The observation of a large 
in-plane resistivity anisotropy, at least for the electron- 
doped 122 Fe arsenides, bears witness to the orthorhom- 
bicity of the material, but does not distinguish between 
anisotropy in the electronic structure and anisotropy in 
the scattering rate. To this end, reflectivity measure- 
ments of detwinned single crystals using polarized light 
can provide important insight to the effects of the mag- 
netic and structural transitions on the anisotropic charge 
dynamics and the electronic band structure. Indeed, 
the counterintuitive anisotropic behavior of p{T) is also 



reflected in the finite frequency response of the charge 
carriers as observed by the optical measurements re- 
ported in our previous Letter work [17 . Optical measure- 
ments of detwinned single crystals of Ba(Fei_a;Co2;)2As2 
in the underdoped regime reveal large changes in the 
low-frequency metallic response on cooling through T^ 
and T/v together with a pronounced optical anisotropy 
(i.e., Aai{uj) = ai{uj,E \\ a) —ai{uj,E \\ b)) at high 
frequencies, defining the linear dichroism [17 . For light 
polarized in the antiferromagnetic a direction, there is an 
increase in the scattering rate, but this is accompanied 
by a dramatic increase in the metallic spectral weight 
that ultimately leads to a reduction in the dc resistivity, 
consistent with observations. For light polarized along 
the b direction, the dominant effect is a reduction in the 
metallic spectral weight, consistent with the increase in 
the dc resistivity. The high frequency dichroism, which 
is smaller for higher Co concentrations, clearly reveals 
that changes in the electronic structure are not confined 
to near the Fermi energy. Similar to dc transport and 
ARPES measurements, a pronounced optical anisotropy 
persists at temperatures above T^ for crystals held under 
uniaxial stress, also anticipating a substantial nematic 
susceptibility. 




FIG. 1: (color online) Sample- holder set-up including a me- 
chanical clamp (a) and an optical mask (b and c). Uniaxial- 
pressure is applied parallel to the 6-axis of the sample (S) by 
drawing the clamp against the side-surface of the crystal (a 
and c). The optical mask, attached in tight contact with the 
clamp device, shapes identically incident and reflected light 
beams for the tungsten reference mirror M (red arrows) and 
S (light blue arrows). 

the full-potential linear augmented plane-wave method 

(LAPW) [m. 



II. EXPERIMENTAL 



In the present work, we expand in greater detail our 
broad spectral range characterization of the anisotropic 
optical conductivity of detwinned Ba(Fei_ajCoa^)2As2 
single crystals under uniaxial pressure. We extensively 
study the underdoped region at compositions with x = 0, 
X = 0.025 and x = 0.045 with light polarized along the in- 
plane orthorhombic a and b axes. We present here the full 
set of data and their complete phenomenological analy- 
sis. Furthermore, in order to clarify the subtle interplay 
of magnetism and Fermi surface topology we compare 
directly our optical measurements with theoretical cal- 
culations obtained from density functional theory within 



A. Samples 

Single crystals of Ba(Fei_a^Coaj)2As2 with x = 0, 2.5 
and 4.5% were grown using a self-flux method [9 . The 
crystals have a plate-like morphology of 0.1^0.3 mm 
thickness, with the c-axis perpendicular to the plane of 
the plates. Crystals were cut in to a square shape, ap- 
proximately 2 mm on the side, oriented such that below 
Ts the orthorhombic a/b axes are parallel to the sides of 
the square [9 . Detailed thermodynamic, transport and 
neutron scattering measurements for the studied dopings 
of Ba(Fei_a;Coa;)2As2 give evidence for structural, mag- 



netic and superconducting phase transitions occurring at 
different temperatures [19l [20]: for x = 0, the coinci- 
dent structural (tetragonal-orthorhombic) and magnetic 
transitions where the system forms antiferromagneticahy 
ordered stripes occur at Tg = T/v = 135 K, whereas for 
X = 0.025 they develop at T^ = 98 K and Tn = 92 K, 
respectively. The compound with x = 0.045 undergoes 
first a structural transition at Tg = 66 K then a magnetic 
transition at Tn = 58 K and finally a superconducting 
one at Tc = 15 K. 

B. Technique 

It has recently been shown that almost single-domain 
specimens can be achieved by application of uniaxial 
pressure in situ [9 . This is crucial in order to reveal 
the intrinsic anisotropy of the orthorhombic phase. To 
this goal we have extended the basic cantilever concept 
originally developed for transport measurements to allow 
optical measurements under constant uniaxial pressure. 
The device consists of a mechanical clamp (Fig. flla)) 
and an optical mask (Fig. [lib) and flic)) attached on 
top of it in tight contact. The pressure-device was de- 
signed according to the following specific criteria: i) the 
uniaxial stress is applied to the sample (S) by tightening a 
screw and drawing the clamp against the side of the crys- 
tal (Fig. flla) and Fig. Illc)). Even though our clamp 
set-up still lacks of a precisely tunable pressure, the uni- 
axial stress was gradually increased, so to observe optical 
anisotropy. The applied pressure is modest, such that T/v 
is unaffected, and can be adjusted over a limited range 
(up to approximately 5 MPa [9 ). Cooling samples under 
such uniaxial stress results in a significantly larger pop- 
ulation of domains for which the shorter ferromagnetic 
6-axis is oriented along the direction of the applied stress, 
almost fully detwinning the crystals, ii) The major axis 
of the tightening screw lies nearby and parallel to the sur- 



face of the sample so that the shear- and thermal-stress 
effects are minimized. The thermal expansion AL of the 
tightening screw, exerting the uniaxial pressure, can be 
estimated to be of the order of AL = aLdT = 20/im (for 
screw-length L = 5 mm, typical metallic thermal expan- 
sion coefficient a = 2 • 10~^K~^, and thermal excursion 
dT = 200 K). This corresponds to a relative variation 
of about 0.4%. By reasonably assuming AL/L = Ap/p, 
the influence of thermal expansion effects is then negli- 
gible, iii) The clamp set-up leaves the (001) facet of the 
single domain samples exposed, enabling us to perform 
optical reflectivity measurements {R{uj)). iv) The optical 
mask guarantees data collection on surfaces of the same 
dimension for S and reference mirror (M) and therefore 
on equivalent flat spots. 

The reflectivity {R{uj)) at room temperature was first 
collected from different spectrometers such as the Bruker 
IFS48 for the mid-infrared (MIR, 500-4000 cm-^) and 
near-infrared (NIR, 4000-7000 cm~^) measurements and 
a PerkinElmer Lambda 950 capable to measure absolute 
reflectivity from NIR up to the ultra-violet (UV) range, 
i.e. 3200-4.8x10^ cm~^. The detwinning device was then 
placed inside our cryostat and finely aligned with micro- 
metric precision within the optical path of the Fourier 
transform infrared Bruker Vertex 80 v interferometer, so 
that we could perform optical measurements of {R{cjo)) 
at different temperatures in the spectral range from the 
far-infrared (FIR, uo < 400 cm"^) up to the MIR, i.e. be- 
tween 30 and 6000 cm~^. Light in all spectrometers was 
polarized along the a and b axes of the detwinned sam- 
ples, thus giving access to the anisotropic optical func- 
tions. 

The real part ai{uj) of the optical conductivity was ob- 
tained via the Kramers-Kronig transformation of R{uj) by 
applying suitable extrapolations at low and high frequen- 
cies. For the uj ^ extrapolation, we made use of the 



3 



1.10 



1.05 



1.00 



0.95 



0.90 



Ba(Fei.xCo, 


)2AS2 


J- ^\ 




\ 








/^^-----:.:.^.-r-50 


^^^:^ 


^-s_- — -^ ^'"^^x-0 


lOK 


x = 


150K 


-■-■ x = 0.025 


lOK 


-■-■ x = 0.025 


150 K 


Cu 

1 1 1 


250 K 

1 



1000 



2000 3000 

Frequency (cm"^) 



4000 



FIG. 2: (color online) Reflectivity ratio R(oj) E\\a I ^(}^^ E\\h 
measured above and below Ts along the a- and 6-axis po- 
larization directions for Ba(Fei-xCox)2As2 with x = and 
0.025 compared with the same ratio for a Cu sample of equiv- 
alent surface dimensions and thickness. 



Hagen-Rubens (HR) formula {R{yj) = 1 — 2. /^), insert- 
ing the dc conductivity values {(Jdc) from Ref. [9 , while 
above the upper frequency limit R{uj) ^ uj~^ {2 ^ s ^ 4) 

Several precautions were taken in order to avoid exper- 
imental artifacts: i) The polarizers chosen for each mea- 
sured frequency range have an extinction ratio greater 
than 200, thus reducing leakages below our 1% error 
limit, ii) As control measurements for the detwinning 
setup we collected at different temperatures the optical 
reflectivity of a Cu sample of comparable surface dimen- 
sions and thickness with respect to the pnictide crystals 
and under equivalent uniaxial pressure. As expected, we 
could not observe any polarization dependence of the Cu 
reflectivity from room temperature down to 10 K (see e.g. 
the data at 250 K in Fig. ^. The Cu test measurements 
set to about 1-2% the higher limit of the polarization de- 
pendence due to any possible experimental artifacts (i.e. 
bended surfaces, leakage of the polarizers etc.), which is 
notably lower then the anisotropy ratio measured for the 
iron-pnictides (Fig. 2 and 3). iii) Prior to performing 



optical experiments as a function of the polarization of 
light, the electrodynamic response of the twinned (i.e., 
unstressed) samples was first checked with unpolarized 
light, consistently recovering the same spectra previously 
presented in Ref. [22]. vi) We achieved the same align- 
ment conditions of M and S (Fig. [lie)) by imaging on 
both spots a red laser point source. 

III. RESULTS 
A. Reflectivity 

The three investigated compositions display overall 
similar features in their optical response but their po- 
larization and temperature dependences show small but 
significant differences as we clarify in the presentation 
and discussion of our results below. Figure |3] presents the 
optical reflectivity R{uj) in the whole measured frequency 
range of detwinned Ba(Fei_a;Coa;)2As2 (x = 0, x = 0.025 
and X = 0.045) at different temperatures and for the 
two polarization directions E||a and E\\b. As already 
recognized in the twinned (i.e., unstressed) specimens 
[22], R{uj) gently increases from the UV to the MIR re- 
gion displaying an overdamped-like behavior. Below the 
MIR energy range, R{uj) gets progressively steeper with 
a sharp upturn at frequencies lower than 200 cm~^ (Fig. 
|3|. Close to zero frequency, R{uj) consistently merges 
in the HR extrapolations calculated with the adc values 
from Ref. [9 . For all measured dopings we observed a 
polarization and temperature dependence of R{uj) from 
the FIR up to the MIR-NIR range, while between 5000 
and 6000 cm~^ the R{uj) spectra tend to merge together. 
The optical anisotropy is rather pronounced at low tem- 
peratures in the FIR region, when approaching the zero 
frequency limit. Interestingly for x = 0, R{uj) increases 
with decreasing temperature along the a-axis, in agree- 
ment with the metallic character of the dc transport prop- 
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FIG. 3: (color online) Temperature dependence of the optical reflectivity of detwinned BaFe2As2, Ba(Feo.975Coo. 025)2 As2 and 
Ba(Feo.955Coo.o45)2As2 in the whole measured spectral range for two polarizations, parallel to the a-axis direction (E||a, red 
solid line) or perpendicular to it (E||6, blue dashed line). The black arrows indicate R(uj) at the temperature close to the 
structural phase transition at T^. The dashed-dotted line (top panels) and the red arrows (bottom panels) indicate the center 
of the MIR band at all temperatures and close to Tn , respectively. 



erties [9 . On the contrary, at low temperatures along the 
6-axis, there is first a depletion of R{uj) in the FIR energy 
range below 700 cm~^ and then a steeper upturn below 
100 cm~^, consistently merging with the HR extrapola- 
tion. While somehow common to all compositions this 
depletion is less pronounced at higher doping levels such 
as X = 0.025 and x = 0.045, when entering the magnetic 
ordered phase at T < T/v (Fig. [sb andlsfc, respectively), 
similarly to what has been previously observed in the 
twinned specimens [22j. 

The anisotropy, discussed so far in the dc-\imit of R{uj) 
at low temperatures, persists up to temperatures well 
above all phase transitions for the crystals held under 
uniaxial pressure, predominantly in the MIR-NIR region 
than in the FIR one. The black arrow in Fig. 3 high- 
light the R{uj) spectra collected at temperatures close 



to the respective T5. The MIR-band centered at about 
1500 cm~^ (dotted-dashed lines and red arrow in Fig. 
3) is of particular interest for both its temperature and 
polarization dependence. For x = 0, we observe an inter- 
change in the polarization dependence of R{uj) with de- 
creasing temperature (Fig. ^ bottom panels). Such an 
interchange exactly occurs at the coupled magnetic and 
structural phase transition at 135 K. At higher doping 
levels {x = 0.025 and x = 0.045) the MIR-band observed 
at 1500 cm~^ for x = is shifted towards lower frequen- 
cies (red arrows in bottom panels of Fig. [3]). Above 150 
K the anisotropy of R{(jj) in the MIR region is strongly 
reduced, but differently from the parent compound the 
a-axis spectrum of x = 0.025 and x = 0.045 is constantly 
above the 6-axis one at all temperatures (Fig. pp andpp 
bottom panels). 
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FIG. 4: (color online) Temperature dependence of the optical conductivity of detwinned BaFe2As2, Ba(Feo.975Coo.o25)2As2 
and Ba(Feo.955Coo. 045)2 As2 in the whole measured spectral range for two polarizations, parallel to the a-axis direction (E||a, 
red solid line) or perpendicular to it (E||6, blue dashed line). The black arrows indicate g\(u) at the temperature close to the 
structural phase transition at Ts , while the red arrows (bottom panels) indicate the center of the MIR-band close to T/v . The 
dashed-dotted and dashed double-dotted line (top panels) mark the frequency cji and uj^ (see text). 



B. Optical Conductivity 

Figure [4] shows the real part G\{ijj) of the optical 
conductivity of detwinned Ba(Fei_a;Coa^)2As2 (x = 0, 
X = 0.025 and x = 0.045) for different measured temper- 
atures along both polarization directions. In the visible 
and UV energy interval (Ti{uj) is characterized by polar- 
ization independent broad absorption bands which over- 
lap with a dominant NIR contribution peaked at about 
5000 cm~^. Consistently with previous measurements 
on twinned samples [22 , these components of ai{uj) are 
generally ascribed to the electronic interband transitions. 
These features in the UV-NIR range are not substantially 
altered by changing temperature, polarization, or doping 
level. 

As already observed for R{uj), the temperature and 
doping dependent optical anisotropy in (Ji{uj) is mainly 



evident in the FIR and MIR regions. In the FIR region, 
there is a strong polarization dependence of the itiner- 
ant charge carriers contribution to (Ji{uj). Along the a- 
axis (Ji{oj) shows a more pronounced metallic behavior 
which gets enhanced below Tn- Along the 6- axis ai{uj) 
below Tn is depleted due to the formation of a pseudo- 
gap, prior to displaying a metallic-like upturn for a; ^ 0. 
As formerly observed for the twinned samples [22 this 
depletion of (Ji{oj) along the 6-axis is less evident when 
increasing the doping. 

The strong absorption peak dominating <Ji (uj) at about 
5000 cm~^ develops into a pronounced shoulder on its 
MIR frequency tail at about 1500 cm~^. As anticipated 
in the presentation of the R{uj) data, this latter MIR- 
band in ai{uj) shows a strong polarization and doping 
dependence, as highlighted in Fig. 4 (bottom panels). 



One recognizes again the already mentioned interchange 
in the polarization dependence when crossing the struc- 
tural transition for x = and its absence at higher dop- 
ings. Above T/v in the MIR range, (Ji{uj) of x = 0.025 
and X = 0.045 along the a-axis direction is above the 
6-axis values and remains constantly above it also below 
T/v (top and bottom panels of Fig. [4t) and [4fc, respec- 
tively). Interestingly enough, for increasing doping the 
maximum of the MIR-band shifts to lower frequencies 
indicating that the MIR-band is significantly affected by 
the doping (dotted-dashed lines in top panels and red 
arrows in the bottom panels of Fig. [4]). 



The anisotropy in the optical response for the magnetic 
state can be anticipated by ab — initio calculations based 
on density-functional-theory (DFT) as well as dynamical 
mean-field theory (DMFT) [HI EH [24! • It was first shown 
that the optical anisotropy of the magnetic state, not 
present within the local spin density approximation, may 
result from DMFT-correlation [23]. Alternatively, DFT- 
calculations of the optical conductivity within the full- 
potential linear augmented plane-wave (LAPW) method 
reproduce most of the observed experimental features, in 
particular an anisotropic magnetic peak located at about 
0.2 eV (1600 cm~^), which was ascribed to antiferromag- 
netically ordered stripes [18 . The optical anisotropy, as 
observed experimentally, was even shown to agree with 
the solution of a three dimensional five-orbital Hubbard 
model using the mean-field approximation in the pres- 
ence of both orbital and magnetic order [25 . Moreover, it 
has been recently pointed out that interband transitions, 
whose relevance is manifested by first-principle calcula- 
tions, give a non negligible contribution already in the 
infrared region, spanning the experimental energy inter- 
val of the MIR-band [26 . We will return later on to the 
comparison between experiment and theory. 
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FIG. 5: (color online) Optical conductivity along the 6- axis 
at 10 K of detwinned Ba(Feo.975Coo.o25)2As2 compared with 
the total Drude-Lorentz fit (black dashed line) and the cor- 
responding components: the narrow (Dn) and broad (D^) 
Drude terms, the mid-infrared (MIR) band and the oscilla- 
tors (/i, /2, h) fitting the interband transitions. The inset 
shows the comparison between the measured reflectivity R{(jj) 
and the resulting flt. 

C. Fits 

In order to study the various contributions shaping 
the optical conductivity at different energies, we apply 
the well-established phenomenological Drude— Lorentz 
approach. Consistent with our previous investigations 
on twinned samples [22 , we ascribe two Drude contribu- 
tions (one narrow and one broad) to the effective metallic 
part of cri(co') and a series of Lorentz harmonic oscillators 
(h.o.) for all excitations (phononic and electronic) at fi- 
nite frequencies. Figure [5] presents all the fitting compo- 
nents for X = 0.025 along the 6-axis direction measured at 
10 K, acting here as a representative example. The use of 
two Drude components in the fit procedure phenomeno- 
logically mimics the multi-band scenario and implies the 
existence of two electronic subsystems as revealed for a 
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FIG. 6: (color online) Temperature and polarization dependence of the widths for the broad (Fb) and narrow (Fat) Drude terms 
(top panels) and the spectral weight (SW) (bottom panels) for the Drude, the MIR and the interband components compared 
with the total SW obtained from the fits of BaFe2As2, Ba(Feo. 975 Coo. 025)2 As2 and Ba(Feo. 955C00. 045)2 As2. The dotted and 
dotted-dashed lines indicate the structural and magnetic transitions at Ts and Tat, respectively. 



wide range of iron-pnictide compounds [27 . The nar- 
row Drude term is relevant at very low frequencies and 
it is obviously tied to the necessary HR extrapolation of 
Riuo) for uo ^0. The broad one acts as a background 
of cFi{uj) and dominates the optical conductivity up to 
the MIR energy interval. As we shall elaborate later on, 
both Drude terms contribute to the total dc conductiv- 
ity. Besides the Drude terms, we chose one broad h.o. for 
the temperature dependent MIR-band and three broad 
h.o.'s (/i, I2 and /a in Fig. Isl) for the strong absorption 
featuring the broad peak centered at about 5000 cm~^. 
The complex dielectric function e = ei{uj) + i£2{oj) can 
be expressed as follows: 
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where £00 is the optical dielectric constant, cjp^, Up^ 
and Fat, F^ are respectively the plasma frequencies, de- 



fined as Up = ^^^^ , and the widths of the narrow and 
broad Drude peaks. The latter parameters represent the 
scattering rates of the itinerant charge carriers, of which 
n, m* and e are then the density, the effective mass 
and the charge, respectively. The parameters of the jth 
Lorenz h.o. as well as those of the MIR-band are: the 
center-peak frequency {uj and uomir)^ the width (7^ and 



Iiir)- The fit 



Imir) and the mode strength (5^ and S 
constraints are such that the measured reflectivity and 
the real part of the optical conductivity are simultane- 
ously reproduced by the identical set of fit-parameters 
[22 , which reduces the degree of freedom in the parame- 
ters choice. The upper boundary for the temperature de- 
pendence of the optical conductivity is found to be close 
to the NIR peak in (Ji{(jj) at about 5000 cm~^. Thus, 
for all temperatures and dopings we fit R{uo) and cri(co') 
by varying the parameters for both Drude terms, the 
MIR-band and the /i h.o., while keeping constant the 
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parameters associated to the two high frequency oscil- 
lators I2 and /a- We systematically adopted our fitting 
procedure (Fig. 5) for both polarization directions and 
for all measured temperatures from 10 K to 270 K of the 
studied compositions. The only exception, however, is 
for the MIR-band of x = along the 6- axis, where we 
added one more component, in order to achieve a better 
fit quality. We checked that this fit- variation has negligi- 
ble impact on the overall trend of the extracted param- 
eters. The remarkable agreement of the fitting results 
with the measured reflectivity and optical conductivity 
(see e.g. Fig. Isl) further demonstrates the overall good 
quality of the fits. Therefore, we are confident that the 
fit results allow us to identify robust trends in relevant 
physical parameters. 

IV. DISCUSSION 

We can first exploit our phenomenological fits from 
the perspective of a so-called spectral-weight {SW) anal- 
ysis, represented by the squared plasma frequencies and 
mode strengths of the Drude terms and Lorentz h.o.'s, re- 
spectively. The Drude-Lorentz procedure thus allows us 
to disentangle the redistribution of SW in selected en- 
ergy intervals and tells us how the same SW is reshuffled 
among the various components as a function of temper- 
ature. Of particular interest is the total Drude weight 
given by SW Drude = ^PN + ^PB' Furthermore, we con- 
sider the SW encountered in the MIR-band and elec- 
tron interband transitions, defined as SWmir = ^'mir 
and SWint = ^j^i^'^i respectively. The total spectral 
weight is then given by the area under the conductiv- 
ity spectrum and can be expressed as [2T1: SWrotai = 
SW Drude + SWmir + SWint- Duc to the high fit-quality 
of cFiiuo) for the chosen components, SWrotai is also 
equivalent to j^"" (Ji{(jj)duj^ where uOc corresponds to a 
cutoff frequency and basically sets the upper frequency 



limit of our measurements. This implies that, if SWrotai 
does not change with temperature within the energy in- 
terval between zero and cjc, any redistribution of the 
spectral weight will fully occur among the fitting com- 
ponents. 

The lower panels of Fig. [6] show the temperature 
dependence of SWrotai and the spectral weight associ- 
ated to the two Drude terms, the MIR-band and the 
three interband excitations (7^). SWrotai stays constant 
with temperature along both polarization directions for 
all doping values thus satisfying the /-sum rule |2T] . 
When crossing the structural and magnetic transitions, 
the spectral weight of the parent compound along the a- 
axis is redistributed from high to low frequencies, since 
the interband transition components loose weight in fa- 
vor of the Drude and MIR-band (bottom panel of Fig. 
|6k). On the contrary along the 6-axis, the Drude term 
looses SW for T < Tat in favor of the MIR-band. At 
higher dopings {x = 0.025 and 0.045) the overah SWjnt 
shows less pronounced temperature dependence for both 
lattice directions. The spectral weight variations with 
temperature mainly occur between the Drude terms and 
the MIR-band for both directions. The resulting tem- 
perature dependence of SWDrude^ increasing along the 
a- but decreasing along the 6-axis across T/v and T5 
for all dopings, is rather intriguing and unanticipated, 
the significance of which with respect to the dc trans- 
port properties will be addressed below. For the time 
being, it is worth emphasizing that the Drude weight 
anisotropy may strongly depend on the topology and 
morphology of the reconstructed Fermi surface below T5 
(i.e., anisotropy of the Fermi velocity), as evinced from 
a five-band Hamiltonian at the mean- field level [28 . Of 
interest is also for x = the temperature independent 
total Drude weight, which is larger along the 6- axis than 
along the a-axis for T > Tg^ thus inverting the polar- 
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FIG. 7: (color online) Temperature dependence of the dichroism Acri(cj) of BaFe2As2, Ba(Feo.975Coo.o25)2As2 and 
Ba(Feo.955Coo.o45)2As2 at uoi and UJ2 (Fig. 4l compared to ^ obtained from the dc transport data, as well as from the 
Drude terms in cfi{uj) {Ap^^^/p ). The vertical dotted-dashed, dotted and double-dotted-dashed lines mark the structural, 
magnetic and superconducting phase transitions at Ts, Tn and Tc, respectively. 



ization dependence observed below Tg. This astonishing 
behavior for x = might be compatible with a recent 
multi-orbital model [25 . Given the ARPES results pT] . 
showing that the band splitting diminishes with increas- 
ing temperature, one would have eventually anticipated 
that the Drude weight gets indeed isotropic above T^. 
Nonetheless, the Drude weight anisotropy above Tg is 
suppressed upon doping (Fig. [6|. Experimentally, such 
a trend of SWorude remains to be verified under con- 
trolled uniaxial pressure conditions, while theoretically 
it awaits confirmation within doping-dependent models. 

From the width at half maximum of the Drude reso- 
nance we can extract the scattering rates of the itinerant 
charge carriers. The broad (T^) and narrow (Tn) Drude 
widths for x = 0, 0.025 and 0.045 are shown in the top 
panels of Fig. [6) An anisotropy in the scattering rate 
is most evident for x = and x = 0.025. For x = 
the Drude scattering rates Fat and F^ increase along 
the a-axis and decrease along the 6-axis for decreasing 
temperatures across the phase transitions. Above Tg for 
X = 0, Fat along both axes saturates to the identical 
constant value, while F^ displays an inversion in the po- 



larization dependence with respect to the situation below 
Tg (Fig. 6, upper panel), before saturating to tempera- 
ture independent values. For x = 0.025 only F^ along 
the a-axis undergoes a sudden incremental change at T/v, 
while all other scattering rates remain almost constant. 
Tb above Tg is switched with respect to below Tg be- 
tween the two polarization directions, similarly to x = 
yet less pronounced. For x = 0.045 the F^ scattering 
rates display a more metallic behavior with decreasing 
temperature above T/v and an almost negligible polar- 
ization dependence. Below Tn there is a weak upturn to 
higher values, particularly along the a-axis. F^r remain 
constant at all temperatures and for both directions. 

The overall temperature dependence of the scattering 
rates below T/v as evinced from the analysis of the optical 
response is expected with respect to the well-established 
magnetic order [4 . Particularly for x = and x = 0.025, 
the larger scattering rates along the elongated antifer- 
romagnetic a-axis than along the shorter ferromagnetic 
6-axis for temperatures below T/v may arise because of 
reduced hopping or of scattering from spin-fluctuations 
with large momentum transfer (i.e., by incoherent spin 
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waves) [29l |30] . The anisotropic scattering rate in the 
paramagnetic state, at least for x = and x = 0.025, 
might also be in agreement with predictions based on in- 
terference between scattering by impurities and by criti- 
cal spin fluctuations in the Ising nematic state [31 . Sim- 
ilarly to our previous discussion on the Drude weight, we 
shall caution the readership, that the trend in the scat- 
tering rates, particularly above Tg, as well as their doping 
dependence should be also verified with tunable uniax- 
ial pressures, in order to guarantee equal experimental 
conditions. A comprehensive theoretical framework, ap- 
proaching different temperature regimes and considering 
the impact of doping, is also desired. 

Having determined the two parameters governing the 
dc transport properties, it is worth pursuing at this point 
the compelling comparison between the temperature de- 
pendence of the optical anisotropy and the anisotropy ra- 
tio of the dc transport properties, defined as —P = Ab^^p^ 
(Fig. ItI) [29 . From the Drude terms, fitting the effec- 
tive metallic contribution of (Ji{uj) over a finite energy 
interval, we can estimate the dc limit of the conductiv- 
ity (a^' = {uj^f/AnTN + {uj^Y/^tiTb) more precisely 
than simply extrapolating (Ji{uo) to zero frequency. The 
anisotropy ratio —^ — , reconstructed from the optical 
data, is thus compared in Fig. [7|to the equivalent quan- 
tity from the transport investigation. The agreement in 
terms of — ^ between the optical and dc investigation is 
outstanding for x=0.025 and 0.045 at all temperatures. 
For x= 0, — ^ — is nonetheless slightly larger than the 
dc transport anisotropy for T < Tg. This disagreement 
might originate from a difference in the applied stress in 
the optical and dc transport measurements, or from dif- 
ferences in scattering rate of samples used for the two 
types of measurements. 

Significantly, analysis of the optical properties for 
all compositions seems to indicate that anisotropy in 



the Fermi surface parameters, such as the enhance- 
ment (depletion) of the total Drude spectral weight oc- 
curring along the a(6)-axis, outweighs the (large at some 
compositions) anisotropy in the scattering rates (Fig. IgI) 
that develops below Tn in terms of the effect on the dc 
transport properties (Fig. [t]). This is an important re- 
sult from the optical investigation, which indeed enables 
to extract both pieces of information governing the be- 
havior of the dc transport properties. 

In order to emphasize the relevant polarization depen- 
dence at high frequencies, we calculate the linear dichro- 
ism Acri(cj), as defined in the Introduction. For the pur- 
pose of further enhancing the optical anisotropy, we show 
in Fig. tel A(Ti(co',T), which is defined as Aai{uj^T) from 
the MIR to the UV for x = 0, 0.025 and 0.045 at various 
temperatures after having subtracted its corresponding 
room temperature values and being appropriately nor- 
malized. The dichroism persists above T/v in the MIR 
range (Fig. Isl), pairing our direct observations in terms 
of R{uj) and ai{uj) (Fig. [3]and|4|. This representation 
highlights once more that the MIR-feature moves towards 
lower frequencies upon increasing doping. 

It is especially interesting to compare the tempera- 
ture dependence of the dc {—^) [9 and optical {Aai{uj)) 
anisotropy [17 . Two characteristic frequencies, iden- 
tifying the position of the peaks in cri{uj) (Fig. [4]), 
are selected to follow the temperature dependence of 
Acri(a;); namely, cji=1500 cm~^ and a;2=4300 cm~^ for 
X = 0; a;i=1320 cm"^ and a;2=5740 cm"^ for x = 0.025; 
co'i=912 cm~^ and a;2=5182 cm~^ for x = 0.045. It is re- 
markable that the temperature dependence of A(Ji{uj) at 
uji and co'2 follows the temperature dependence of — ^ in 
both compounds (Fig. [7|. Acri(a;^) (z = 1, 2) saturates at 
constant values well above Tg and then displays a varia- 
tion for T < 2T5. Here we first underscore that the rather 
pronounced optical anisotropy, extending up to temper- 
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FIG. 8: (color online) Temperature dependence of the normalized relative difference Aai (cj, T) = Aai (cj, T) — Aai (cj, 300 K) of 
the optical dichroism of BaFe2As2, Ba(Feo.975Coo.o25)2As2 and Ba(Feo.955Coo.o45)2As2. The black arrows indicate the Aai{(jj) 
at the temperature close to the structural phase transition at T^. As color code: red marks positive, green negative values of 

Acri(cj,T). 



atures higher than Tg for the stressed crystals, clearly 
implies an important pressure-induced anisotropy in the 
electronic structure, which is also revealed by ARPES 
measurements [61ITT]. 

Since the dichroism directly relates to a reshuffling of 
spectral weight in ai{uj) in the MIR-NIR range (Fig. 4a 
and 4b), A(Ji{uj) at uji is interrelated to that at UJ2 (i.e., 
the right ?/-axis for Aai{uJi) in Fig. [tI are inverted be- 
tween ui and CJ2), so that the behavior of Aai{uj) is 
monotonic as a function of temperature and opposite in 
sign between uji and 002 (Fig. [7|and|8|. For x=0.025 (Fig. 
[7| Aai{uji)=0 at T >> T^. However for x= and 0.045, 
A(Ji{uJi) is found to be constant but apparently different 
from zero for T > > Tg . The origin of the finite (but con- 
stant) dichroism at high temperatures for these samples 
is at present unclear, and might reflect a systematic effect 
due to imperfect experimental conditions (e.g., too strong 
applied uniaxial pressure). Nevertheless, the overall tem- 
perature dependence seems to behave in a very similar 
manner for all compositions. Significantly, the absolute 



variation of the dichroism across the transitions at se- 
lected frequencies is larger for x=0 than for x=0.025 and 
0.045, contrary to the anisotropy in the dc resistivity [9]. 
This doping-dependence needs to be studied in a con- 
trolled pressure regime in order to exclude effects arising 
from different degrees of detwinning (T < Tg) and dif- 
ferent magnitude of induced anisotropy (T > Tg). Even 
so, it is encouraging that, contrary to the dc resistivity, 
the changes in the electronic structure appear to follow 
a similar trend to doping as the lattice orthorhombicity 
[T5]. Our data might thus reveal a pronounced sensitiv- 
ity of the electronic properties to structural parameters, 
like the iron-pnictogen angle a [32 , altered by external 
tunable variables like uniaxial pressure. Indeed, changes 
in a seem to induce relevant modifications in the shape 
of the Fermi surface and its nesting properties as well 
as in its orbital makeup, thus implying consequences in 
terms of the superconducting order parameter, critical 
temperature and magnetic properties ^32] . 

The origin of the orthorhombic transition has been dis- 
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cussed from the closely related perspectives of spin fluc- 
tuations (a so called spin-induced nematic picture [33]- 
[36]), and also in terms of a more direct electronic ef- 
fect involving, for instance, the orbital degree of freedom 
[25[l28ll29l[37H4Q] . The present measurements alone can- 
not distinguish between these related scenarios, because 
in all cases some degree of electronic anisotropy is an- 
ticipated, and indeed it is likely that both orbital and 
spin degrees of freedom play a combined role in the real 
material. Nevertheless, it is instructive to compare the 
observed dichroism with specific predictions made within 
models based on orbital order. The two well-defined en- 
ergy scales uji and cc;2 (Fig. [4| may represent optical tran- 
sitions between states with the strongest d^zl^yz charac- 
ter, which are separated by about 0.3-0.4 eV. Such an 
energy splitting is indeed compatible with the theoret- 
ical calculations of the anisotropic optical conductivity 
[25] and of the linear dichroism in the X-ray absorption 
spectroscopy [29 . Nonetheless, the debate about the im- 
pact of orbital order on the electronic properties is far 
from being settle down. Valenzuela et al. point out 
that an increasing degree of the orbital order may fa- 
vor a diminishing anisotropy of the excitation spectrum 
and above all of the Drude weights [28] , opposite to what 
we may conclude from our optical experiment. Bascones 
et al. even claim that the orbital order accompanies the 
magnetization within a wide range of parameters but it 
is not correlated with the magnetic exchange anisotropy 
[40] . While these latter calculations might be relevant 
at low temperatures, it remains to be seen how they can 
explain the onset of the anisotropy at/above T^, where 
the spin symmetry is not yet broken. 

Finally, we wish to come back to the comparison be- 
tween our optical results and the outcome of the LAPW 
calculations [I8l [41]. The model uses as initial lattice 
parameters those determined experimentally from the 



crystal structure of BaFe2As2 but neglects the small 
orthorhombic distortion occurring at low temperatures. 
The magnetic ordering is modeled using a fixed striped 
phase, where stripes of spins up are alternated with 
stripes of spins down on each Fe plane, compatible with 
the magnetic configuration determined experimentally 
below T/v [4]. The Co doping of BaFe2As2 rapidly desta- 
bilizes the magnetic phase which disappears above 6% Co 
concentrations. This Co-induced destabilization of mag- 
netism was simulated by scaling the magnetic moment, 
calculated within the virtual crystal approximation [18], 
accordingly to the experimental critical temperature T^ 
[To] , This approach corrects the deficiencies of the lo- 
cal spin density approximation and mimics the two main 
Co-induced effects: a localized electron doping on the Fe 
layer and a magnetic to nonmagnetic phase transition. 
The resulting theoretical Co-dopings agree within 0.5% 
variation with the experimental ones. 

Figure [9] show the low temperature measured (top pan- 
els) and calculated (bottom panels) optical conductiv- 
ity for both the a- and 6-axis and for the three dopings. 
For a direct comparison we have normalized all the mea- 
sured and calculated g\(uj) to their respective maxima, 
thus obtaining G\(bj). We clearly see a fairly good agree- 
ment between theory and experiment in the general shape 
of G\(ijj). In the MIR range the DFT calculated ai(co') 
finely reproduces the observed polarization dependence, 
both as far as the MIR-band position and its polariza- 
tion dependence (black arrows in Fig. [9| are concerned. 
Indeed, the predictions of the enhancement of the MIR- 
band along the a-direction and its depletion along the h- 
direction are fairly close the experimental findings. Inter- 
estingly, the center of the calculated band (black arrows 
in Fig. [9] bottom panels) shifts at lower frequencies for 
increasing dopings which is also in agreement with the 
experiments. This MIR-band is linked to the modeled 
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FIG. 9: (color online) Normalized optical conductivity d-i{(jj) of BaFe2As2 (a,b), Ba(Feo.975Coo. 025)2 As2 (c,d), and 
Ba(Feo.955Coo.o45)2As2 (e,f) measured (top panels) at the lowest temperatures with light polarized along the a- and 6- axis 
and obtained by DFT calculations (bottom panels) along the antiferromagnetically oriented spins and parallel to the ferromag- 
netic stripes ^8j. The black arrows indicate the MIR-band in the measurements and calculations. 



magnetic stripe configuration which was shown to cor- 
respond to the energy- minimum configuration of these 
systems [18 . Therefore, the DFT calculation strongly 
supports a "magnetic origin" of the MIR-band which 
would originate from the Fermi topology reconstruction 
in the magnetically ordered state. In this scenario one 
would reasonably expect that the MIR-band disappears 
above the magnetic phase transition, contrary to our ob- 
servations. However, a dynamic antiferromagnetic order 
due to spin fluctuations could persist in the paramag- 
netic phase well above the phase transition temperature 
[42j . The fingerprints of such underlying spin fluctuations 
would be frozen to fast enough probes as optics, thus 
explaining the persistence above the phase transition of 
the MIR-band in our spectra. At higher (NIR) frequen- 
cies the agreement result depleted because of the finite 
/c-point sampling. The major interband peak of d-i{uj) ex- 
perimentally observed at about 5000 cm~^ is shifted to 
slightly higher frequencies in the theoretical calculations 
and shows a steeper rise. We notice that spin-polarized 
DFT calculations require a smaller renormalization fac- 



tor with respect to unpolarized ones in order to account 
for these frequency- shifts. 



V. CONCLUSIONS 

The charge dynamics of detwinned Ba(Fei_ajCoa;)2As2 
single crystals in the underdoped regime reveals an 
in-plane temperature and doping dependent optical 
anisotropy. At low frequencies the optical measurements 
offer the unique opportunity to disentangle the distinct 
behaviors of the Drude weights and scattering rates of 
the itinerant charge carriers, which are both enhanced 
along the antiferromagnetic a-axis with respect to the 
ferromagnetic 6- axis. Our findings on such single do- 
main specimens allow us to shed light on the counterintu- 
itive anisotropic behavior {p^ > pa) of the dc resistivity. 
The dc anisotropy below T/v is principally determined by 
the anisotropy in the low frequency Drude weight (i.e., 
changes in the electronic structure close to the Fermi en- 
ergy), outweighing the non-negligible anisotropy of the 
scattering rates between the a- and 6- axis. Of equal or 
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perhaps greater interest is the temperature regime above 
T/v for which the Fermi surface is not reconstructed. 
One would hke to understand whether the resistivity 
anisotropy at high temperatures also originates from the 
Fermi surface, perhaps due to the difference in the or- 
bital occupancy revealed by ARPES [25l [25, or from 
anisotropic scattering, perhaps associated with incipient 
spin fluctuations [31 . The current optical data may be 
in partial agreement with both points of view and there- 
fore do not permit a conclusive answer to this question, 
thus motivating further experiments in order to definitely 
address the origin of the anisotropy above T/v . 

The optical anisotropy extends to relatively high fre- 
quencies and temperatures above the phase transitions 
for crystals held under uniaxial stress. The resulting 
linear dichroism reveals the electronic nature of the 
structural transition and implies a substantial nematic 
susceptibility. In order to clarify the subtle interplay 
of magnetism and Fermi surface topology we elaborate 
on a comparison of our optical measurements with 
theoretical calculations obtained from density functional 
theory within the full-potential LAPW method. The 
calculations are able to reproduce most of the observed 



experimental features, in particular, to identify the 
MIR-band located at about 1500 cm~^ as a magnetic 
peak, ascribed to antiferromagnetic ordered stripes. 
The measured large in-plane anisotropy of the optical 
response and its doping dependence is consistently 
tracked by the LAPW calculations. 
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